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(Manuscript received November 9, 1956) 

The number of calls in pi-ogress in a sijnple telephone exchange model 
characterized by unlimited call capacity, a general probability density of 
holding-time, and randomly arriving calls is defined as N{t). A formida, 
due to Riordan, for the generating function of the transition probabilities 
of N{1) is proved. From this generating function, expressions for the co- 
variance function of N{t) and for the spectral density ofN{t) arc determined. 
It is noted that the distributions of N{t) are completely specified by the co- 
variance function. 

I INTRODUCTION 

The aim of this papei' is to study the average fluctuations of telephone 
traffic in an exchange, by means of a simple mathematical model to 
which we apply concepts used iu the theory of stochastic processes and 
in the analysis of noise. 

The mathematical model we use is based on the following assumptions: 
(1) requests for telephone service arise individually and collectively at 
random at an average rate of a per second; (2) the holding-times of 
calls are mutually independent random variables having the common 
probability density function h{u)] and (3) the capacity of the exchange 
is effectively unlimited, and no call is blocked or delayed by lack of 
equipment. This telephone exchange model has been described by J. 
Riordan/' 

As a measure of trafhc, it is natural to use the number of calls in prog- 
ress in the exchange. We are thus led to consider a random step-function 
of time N(t), defined as the number of calls in progress at time t. N{t) 
fluctuates about an average in a manner depending on the calling-rate, 
a, and the holding-time density, h{u). 

II PROOF OF RIORDAN '« FORMULA FOR TRANSITION PROBABILITIES 

Let Pm.,Xl) be the probability that n calls arc in progress at ( if m 
calls were in progress atO. Define the generating function of these prob- 
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abilities as 

and let 

fill) = ( h{x) dx, 

SO that the average holding-time, h, is given by 

h = I fi:u) du. 
Jo 

Iliordan' has given the following I'ormula for P^t, x) : 

PUt, 0-) = [1 + (x - l)g(t)r exp {{x - l)ah[l - git)]], (1) 

with 

git) -l[ fiu) du. 

For exponential holding-time density, this formula had ah-eady been 
derived (as the solution of a differential equation) by Palm.' 

In private communication, J. Riordan has suggested that his proof of 
(1) is incomplete. We therefore give a new proof of (1). 

We seek the generating function of Nit), conditional on the event 
]V"(0) = m. We obtain it by first computing the joint generating function 
of iV(0) and N{t); that is, 

£;[/'V*"}- (2) 

The desired conditional generating function is then the coefficient of 
y'" in (2), divided by the probability that iV"(0) = m. 

To obtain a formula for (2), we exhaust the interval (— «, 0) by 
division into a countable set of disjoint intervals, /„ , the h^'' having 
length Tn > 0. Let S„ be the sum of the first n lengths, Tj . Let ^„(/), 
for t > — Sn-i , be the number of those calls which arrive in /„ and are 
still in progress at t. And let rjit) be the number of calls arriving during 
(0, (), t > 0, and still in existence at t. Then 

A^(0) - E ^-.(0), (3) 

Nit) = vit) + E Ut), t > 0. (4) 

Since calls arriving during disjoint intervals are independent, we know 
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that 7}{t) is independent of all the ^'s, and that ^„{i) is independent of 
$j(r) if n 7^ j. Of course, f„{0 and ^„{t) are not independent. It follows 
that if the infinite product converges, then for i > 

/i:{/r-r"'"l =ii^{x"") n£|?/"'V"'"}. (5) 

We now compute the terms of the product. If a call originates in in- 
terval /„ , it still exists at with probability 

I ,, Jo in J8„-i 

Hence if k calls arrived in /„ , the probability that m of them are still 
in progress at is 

prUn(O) = 'm\k calls arrive in /„) 

= (^)q"(1 -Q.Y-"\ m^k. 

Similarly, if a call originates in /„ and exists at 0, it also exists at f > 
with probability 

Kn = {QnTnV [ "/(» + t + S„.i) du. 
Jq 

Therefore 

/i'{a;*"''^ ! $n(0) = m and k calls arrive in /„} 

= [1 + (x - 1)KX> 
and so 

= (1 + (^[1 + {x - 1)K„] - 1)Q„)'- 



k 
= a . 



The number of calls arriving during /„ has a Poisson distribution with 
mean aTn ; hence 

£,{^U(0)^t.cOj =exp{ar„Ca- 1)1 

(6) 

= exp {aTnQ.{y\\ + {x - l)K.] - 1)J. 

By reasoning like that leading to (6), it can be shown that 
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£(.1:"'") = exp Lt{x - 1) ^ ^ /(") du 



<7) 



= exp {ah(x - m - giOW- 



Now 



23 al\Q„ = a I f{u) du = ah, 

= a I /(w) rf« = ahg{i). 

Therefore the infinite product is convergent, and 

E{y'''''x'""\ = exp {ahi:x. - 1)[1 - (j{t)\ 

+ E ar„Q„(i/[l + {x ~ \)K^ - \)] (8) 

= exp Saft.((:r - 1)[1 - s(01 + 0/ - 1) + ?/(-^^ " Df/W)}- 

Thus the generating function of the joint distribution of A^(0) and N{t) 
is independent of the division of {- ^ , 0) into intervals In ■ By letting 
X approach 1 in (8) and finding the coefficient of if in the resulting limit, 
we find that 

pr{iV(0) =m\ ='-^^^- (9) 

The coefficient of -if in (8) itself is 
e-''\ah) 



m\ 



[1 + Or - l)p(«)l"'exp \{x - DaJiW - (/(/)]}, 



and so using (9) we find that the required conditional generating function 
of N{t), given iV(0) = m, is given by Riordan's formula (1). 



Ill THE AUTOCORRELATION 



In terms of N{t) one can define various stochastic integrals which 
will be characteristic of the process. A simple one which has been ex- 
tensively treated in connection with estimating the average traffic is 
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the average of N{t) over an interval (0, T). The chief references in the 
literature on M are References 3 and 5. If we consider N{t) during an 
interval (0, 2' -f t), a measure of the coherence of N{t) during this in- 
terval, i.e., of the extent to which Nit) hangs together, is given by the 
integral 



Ua\ r) ^1= [ NO) N(t + r) dt, 
1 Jo 



depending on values of N(t) taken r apart. When the limit i^(r) of u as 
T appi'oaches ^ exists, it is usually called the autocorrelation function; 
most statisticians, however, reserve the term "correlation" for suitably 
normalized, dimensionless quantities. It can be .shown that this limit 
exists and is the same for almost all A'^(/) in the ensemble. It then coin- 
cides with the ensemble average, i.e., 

i^(r) = lim U{T, t), almost all Nit), 

= E{Nit)Nit-^ t)\. 

The function, ^, for the system we are discussing is derived by Riordan," 
and we reproduce his argument for ease of understanding. For equilib- 
rium, and b = ah, we have 

E{Nit)Nit-\-r)} = j:i^m^P„Xr,x)] . 
m-o ml ax Ji-i 



Now 



so that 



d.r 



r=l 



= mgir) + 6[1 - gir)]. 



Ht) = E^-^m{mgiT) + Ml - g(r)]\, , ^ 

«^ ml (10) 

= h' + hgir). 

(Cf-.-'p. 1136) 

The limiting value of i/'(t) for t approaching x i.s the square of the 
mean occupancy, 6, and the limiting value of ^(t) for t approaching 
is the mean square occupancy, b' + b, the second moment of the Pois- 
son distribution with mean b. 
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IV THE COVAHL\NCE AM) SPECTRAL DENSITY 

The average value of jV(/) isb ^ ah. One way to study the fluctuations 
of Nii) about its average is by means of the power spectrum used in the 
analysis of noise. (Cf. Rice/) We resolve the difference [N{t) - b] into 
sinusoidal components of non-negative frequency, and postulate a noise 
current proportional to this difference dissipating power through a unit 
resistance. The spectrum w{f) is then the average power due to frequen- 
cies in the interval (/, / + df) . 

More formally, we consider the Fourier integral 



Sif, T) = f [Nit) - b]e-'"'' dl, 
Jo 



and we recall, for completeness, the relationship between S{f, T) and 
the covariance function, /?(t), of [N{f) — b]. If 

w(f) =lim I y ' -^, 
then 

- w{f) = 4 f R(t) cos27r/Trfr, 

t (11) 

Rir) = f w(f) cos 2ir.fT df. ' 
Jo 

(d. Rice,' p. 312 ff.) 
At the same time, we have 

R{r) = E{[N{t) - b][N{t-\- r) - b]\ 

= ^(r) - b' 

= bgir). 

Let X{t) be any stochastic process which is known to be the occupancy 
of a telephone exchange of unlimited capacity, having a probability 
density of holding-time, and subject to Poisson traffic. From the pre- 
ceding result it can be seen that the covariance function of X{1) deter- 
mines the distributions of the X{t) process completely, since 

dH 



fir) = [h(u)du = -a''^. 



If the holding-times are bounded by a constant, k, then readings of 
A'^(() taken further apart than k are uncorrected . In fact, such values 
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are independent, because no call which contributes to N{t) can survive 
until (t + k), with probability 1. 
Using (11), we see that 



w 



(/) = 4 f cos 2TrfTR(T) dr 

Jo 

= 46 / cos ^-KfTgir) dr 
Jo 

= 4a / cos 27r/T / / h{u) du dy dr (12) 

= —f I sin 2TrfT / k{u) du dr 
vj Jo Jt 

= ^[l - { C0s2TM(T)(?rl. 

Equation (12) expresses the mean square of the frequency spectrum of 
the fluctuations of the traffic away from the average in terms of the call- 
ing-rate and the cosine transform of the holding-time density, h{u). 
The calling-rate appears only as a factor, and so does not affect the shape 
of w(J). The function «?(/) is what Doob^ (p. 522) calls the "spectral 
density function (real form)." 

V EXAMPLE 1. N(t) MARKOVIAN 

Let the frequency h(u) be negative exponential, so that 

h(u) = I e-"'\ (13) 

where h is the mean holding-time. It is shown in Riordan** p. 1134, 
that N{t) is Marlcovian if and only if h{u) has the form (13). From page 
523 of Doob^ we know that the covariance function of a real, stationary 
Markov process (wide sense) has the form 

Rir) = R{0)e~''\ a constant. (14) 

Under the assumption (13), the covariance of N{t) is 

R{t) = bgir) =1} f h{u) du dy 



= be 



-T/ft 
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in agreement with (14). The spectral density can now be obtained from 
(U) or (12); it is 



w{f) = 



4fa/t 



1 + AwTh^' 



This is the same as would be obtained for a Markov process that alter- 
nately assumed the values -i-y/ah, -\/ah at the Poisson rate of (2/i)"^ 
changes of sign per sec. (Cf. Rice p. 325.) 

VI EXAMPLE 2. HOLDING-TIME DISTRIBUTED UNIFORMLY IN {a,0) 

Let hiu) be constantly equal to (|3 - a)~^ in the interval (a, /3), and 
constantly elsewhere. Then by (12), 

_ a r. sin 27r/^ — sin 2irfa 



Now we see that 



I(y) = f h{u) du = 
•ill 



2irf{0 - a) 

1 for y ^ ct 

/3 - 2/ 



— a 
P for y ^ ^ 



for a S y ^ 



so that 



a |^« - r -H ^^J 



Rir) = 



a {0- rf 
2 - oi 





< T ^ 



a ^ T £ /3 



(15) 



r S 



is the covariance function of the process N{t) when holding-time is dis- 
tributed uniformly in (a, /3). 

If, formally, we let (/3 - a) approach while keeping |(a + /3) fixed, 
then the holding-times become concentrated in the neighborhood of the 
mean, h; in the hmit, as /i(w) tends to a singular normal distribution 
with mean, h, and variance zero, we obtain 






(16) 
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as the spectral density function for the N{t) process with constant hold- 
ing-time, h = ^{a -\- /3). Similarly, from (15), we note that as the hold- 
ing-times become singularly normal -with mean, h, and variance zero, 
the covariance function becomes 

i= aQi - t) ^ t :£ h 

1=0 T > h. 



We can express (16) as 



and note that this is exactly like the power spectrum of a random tele- 
graph wave constructed by choosing values + y/ah, — y/ak with equal 
probability and independently for each interval of length, h. (Cf. Rice,* 
page 327.) 
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